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Abstract 

Cross sections for the disintegration of the deuteron via neutral-current (NCD) 
and charged-current (CCD) interactions with reactor antineutrinos (u e d — > u e pn and 
v e d -> e + nn) are measured to be 6.08 ± 0.77 x 10 -45 cm 2 and 9.83 ± 2.04 x 10 -45 cm 2 
per neutrino, respectively, in excellent agreement with current calculations. Since the 
experimental NCD value depends upon the CCD value, if we use the theoretical value 
for the CCD reaction, we obtain the improved value of 5.98 ±0.54 x 10 -45 for the NCD 
cross section. 

The neutral-current reaction allows a unique measurement of the isovector-axial 
vector coupling constant in the hadronic weak interaction, j3. In the standard model, 
this constant is predicted to be exactly 1, independent of the Weinberg angle. We 
measure a value of j3 2 = 1.01±0.16. Using the above improved value for the NCD cross 
section, /3 2 becomes 0.99 ± 0.10. 



1 Introduction 

We describe an experiment to measure the cross sections for the disintegration of deuterons 
by neutral- and charged-current interactions with low energy electron-antineutrinos. Data 
were taken at the Centrale Nucleaire de Bugey in France, at 18 m from the core of Reactor 5. 

Improvements were made to the cosmic-ray shielding of the detector which we previously 
used in a similar experiment at the Savannah River Plant in South Carolina in the late 
1970s An outer layer of active cosmic ray veto detectors was added which completely 

surrounds the lead and steel gamma ray shield. These improvements reduced the neutron 
background due to cosmic rays by a factor of six to ~ 25 day -1 . 
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There are two reactions of interest in this experiment — the Neutral Current disintegra- 
tion of the Deuteron (NCD), 

v e + d — > v e + p + n, 
and the Charged Current disintegration of the Deuteron (CCD), 

v e + d — > e + + n + n. 

The experiment was designed to probe these reactions at low energies (~1 MeV). In 
particular, it measures the square of the isovector-axial vector coupling constant {(3 2 ). The 
neutrino-induced disintegration of the deuteron is an ideal reaction for this purpose since, 
at reactor neutrino energies, all other coupling constants make negligible contributions to 
the cross section. Other coupling constants depend on the value of the Weinberg mixing 
angle, 9w, which is an unspecified parameter of the theory, while (3 is predicted to be -1.0, 
independent of 6w- In addition, it does not suffer from ambiguities arising from the presence 
of vector interactions, nor from momentum-transfer-dependent form factors, to which high- 
energy experiments are subject. The deuteron disintegration experiment is unique, then, in 
being able to measure the contribution of a single coupling constant with an unambiguous 
theoretical value. 

2 The Detector 

2.1 Location 

The detector was installed at Reactor 5 of the Centrale Nucleaire de Bugey, near Lyon, 
France. It is located in a room about 10 m below ground level with an overburden of 
25 mwe. The distance from the center of the reactor core to the center of the detector was 
18.5 meters. 

2.2 The Target 

Schematics of the detector and shielding are shown in Figures [I], 0, and |3|. The target 
detector, labeled D 2 in Figures [3] and @ and shown in more detail in Figure ||, consists of 
a cylindrical stainless steel tank, 54 cm in diameter, 122 cm in height, and a wall thickness 
of 0.18 cm, containing 267 kg of 99.85% pure D 2 and ten tubular proportional chambers, 
equally spaced in two concentric rings of 10.16 cm and 20.37 cm radius and offset from each 
other by 36°. Running down the center of the target tank is a stainless steel tube that allows 
the placement of radioactive sources inside the tank for calibration purposes. 

Immediately surrounding the target tank is 10 cm of lead shielding and a 1 mm layer of 
cadmium to absorb thermal neutrons. These are contained in an outer steel tank that sits 
on a small pedestal inside the large, inner veto detector tank (Tank 2 of Figures [1] and |2]). 

The proportional counters are 5.08 cm in diameter, 122 cm in height, have a wall thickness 
of 0.025 cm, and are filled with 1 atm of 3 He and 1.7 atm of Ar as a buffer. They are essentially 
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Figure 1: Plan view of detector and shielding in the opened configuration at the Bugey site. 
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Figure 2: Side view of detector and shielding configuration. 




Figure 3: Top view of target tank. 
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Figure 4: Neutron-response spectrum of a 3 He counter with a 252 Cf source. 

black to thermal neutrons, with a capture cross section of ~5300 barns per 3 He nucleus. The 
neutron capture in the counters proceeds via the (n,p) reaction: 

3 He + n -> 3 H + p + 764 keV. 

The energy resolution of the counters was measured to be 3% at the 764 KeV neutron capture 
peak. A typical neutron spectrum obtained with a 252 Cf neutron source is shown in Fig. [|. 
A discussion of the neutron detection efficiency is given in Section |3.3| . A more detailed 
description of the construction and testing of the 3 He proportional tubes can be found in 
Reference |§. 

2.3 Detection Technique 

The neutral-current and charged-current events in the D2O target are recognized solely by 
the neutrons they produce: the neutral-current reaction releases a single neutron and the 
charged-current releases two. Consequently, the quantities of interest are the rates of single 
and double neutron captures. 

2.4 The Shielding and Anticoincidence System 

Due to the detector's close proximity to the reactor core, there can be a significant reactor 
associated gamma flux. Gamma rays of >2.2 MeV which reach the target detector can 
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photodisintegrate the deuterons, leading to single neutron signals. In the previous version 
of the experiment, this background was reduced by surrounding the inner layer of active 
cosmic-ray veto detectors with a layer of lead and water shielding. Unfortunately, cosmic 
rays interacting in the surrounding lead shield, but not reaching the inner veto counters, 
were a significant source of neutrons in the target detector. 

It was concluded that the shielding could be improved by an additional layer of active 
cosmic ray veto detectors outside the lead shielding. In this way, cosmic rays interacting 
in the lead would be seen by the outer veto detectors. Simulations showed that this would 
reduce the cosmic ray neutron background by a factor of three to four. 

In the current configuration, the target tank is in the center of a large liquid scintillator 
detector (the "inner" veto) composed of Tank 2 and Tank 4, shown in Figure |2|. Immediately 
surrounding Tanks 2 and 4 is a layer of each lead and steel. Surrounding this layer of passive 
shielding is an outer layer of cosmic ray veto detectors (the "outer" veto). Slabs of plastic 
scintillator cover the north and south sides and the bottom face, while larger tanks of mineral 
oil scintillator (Tanks 1, 3, & 5) cover the east, west, and top faces. The liquid scintillator 
used in all five tanks is mineral oil based with a high flash point. Five-inch hemispherical 
photomultiplier tubes (PMTs) are used to view the liquid scintillator tanks and three inch 
tubes are employed on the plastic slabs. 

As noted above, the inner veto system consists of two liquid scintillator tanks, Tank 2 
and Tank 4. As indicated in Figures p] and |2j, there is a string of fifteen evenly spaced 
PMTs along each vertical corner of Tank 2 . Alternating tubes are offset in direction by 
90°. Along the east and west walls on the floor of the tank is a row of PMTs which view the 
space underneath the target tank. Tank 4 has three PMTs in each vertical corner that are 
configured like those of Tank 2. 

There are eight signal lines coming from the inner veto system. The PMTs in each 
vertical string are ganged onto a single line, as are each row along the floor. The signals 
from the east corners of Tank 4 are fanned together as are the signals from the west corners. 

The inner veto detector is primarily a "soft" veto — its signals are recorded at each 
trigger and analyzed off line. However, it also triggers an on-line veto in the event that all 
four corner strings see a large pulse simultaneously. Such a signal is likely to be produced 
by a throughgoing muon. 

Figure |^ shows some details of our electronics configuration. 

2.5 Data collection system 

The data-collection program, based on a 80486DX processor and the software package Lab- 
VIEW, has a fast, graphical interface to the electronics. The software takes advantage of 
the multitasking capabilities of the operating system, allowing the transfer and processing 
of data without interrupting data collection. 

A trigger is generated under the following conditions: 

1. A neutron-like pulse is detected in one of the 3 He proportional counters. 
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Figure 5: Some details of our electronics configuration. 
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Figure |] continued. 
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2. No pulses above hardware thresholds were detected in any of the inner or outer veto 
detectors in the preceeding ~ 900/is. This value was chosen to reduce background 
from muon-induced neutrons arising in the inner-anti scintillator which had a neutron 
capture time of about 200/is. 

When a trigger occurs, the contents of waveform digitizers and scalers are read by the 
computer and written to disk. The contents of the digitizers give a pulse history of all 
detectors for a period of 4 milliseconds before and after the event. 

More details on the detector and data-collection system can be found in Reference ||. 

3 Data Analysis 
3.1 Selection criteria 

After the data are collected, they are further reduced by offline selection according to the 
following criteria. 

3.1.1 Target cuts 

The purpose of the target cuts is to remove any events that do not appear to be valid neutron 
captures. 

No neutron in pulse-height window. An event is removed if there is no target pulse in the 
pulse- height acceptance window within 5/isec of the trigger time. The pulse- height 
acceptance values were determined from the neutron calibrations and varied slightly 
from run to run. The total number of target peaks in the pulse-height window during 
the 782/isec (three times the neutron capture time in the target) following the trigger 
is taken to be the number of neutrons in the event. This time interval was selected 
to maximize the signal to background. Shorter time windows yield consistent results 
with larger statistical errors. 

Remove "early" neutrons. Remove event if it has a pulse in the neutron pulse- height accep- 
tance window before the trigger time. 

3.1.2 Outer-anti cut 

This cut removes cosmic-ray muons that might create neutrons that would subsequently be 
detected by the target. 

Remove muons. If during the 1800 /isec preceeding the trigger a signal is detected in any 
outer anti which exceeds the threshold for that counter, the event is removed. 

The pulse-height threshold values were determined run by run. The values were chosen 
at the lowest value for which the events removed by this criterion were at least twice 
the number of the "background" peaks at the same pulse height. 
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Reactor ON 


Reactor OFF 


Cut 


Random 


Normal 


Random 


Normal 




Triggers 


Triggers 


Triggers 


Triggers 


Early neutron 


0.5% 


27.9% 


0.5% 


28.7% 


No neutron in window 


0% 


17.2% 


0% 


17.5% 


Outer-anti cut 


26.0% 


38.9% 


26.4% 


38.7% 


Inner-anti (high-energy) 


9.2% 


8.5% 


9.5% 


9.3% 


Inner-anti (low-energy) 


17.2% 


6.5% 


19.2% 


5.1% 



Table 1: Fractions of total numbers of all events removed by each data cut when applied 
in the order indicated. Since many events satisfy more than one cut criterion, these values 
would change if the ordering of the cuts were changed. 

3.1.3 Inner-anti cuts 

The inner-antis provide additional protection against cosmic-ray muons that sneak through 
the outer anti. However, this large volume of liquid scintillator also provides a large target 
for inverse-beta events on hydrogen (u e p — > e + n). A small fraction (~ 0.08%) of the neutrons 
thus produced diffuse into the target area and are recorded by the 3 He tubes. To reduce this 
number, a low-energy cut is applied to the inner antis, thus using the light produced by the 
positron to veto the event. 

Low-energy cut. The main purpose of this cut is to remove the inverse-beta events. Any 
event with >0.8 MeV in either Tank 2 or Tank 4 within 900/isec before to 200/isec 
after the trigger was removed. This energy-threshold value was chosen in order to 
remove the maximum number of inverse-beta events, while not suffering too much 
dead time from the many low-energy background pulses. From the Monte Carlo, the 
mean time between production of a neutron by the inverse-beta process in Tanks 2 or 
4 and its subsequent capture by a 3 He tube was about 230/isec. Thus the period in 
which this cut is active is from about four such mean times before to one after the 
trigger. 

High- energy cut. Events having a pulse of total energy exceeding 8 MeV in Tank 2 or 6 MeV 
in Tank 4 from 2400 to 900/isec before the trigger are removed. This cut removes 
cosmic-ray events that are recorded before the beginning of the hardware anti block. 
Extending the times earlier than 2400/zsec has little effect. 

The fraction of events removed by each of the above cuts is shown in Table [l]. Figure ^| 
shows the effects of the cuts on the data. As a result of the cuts, the number of candidate 
neutrino events is reduced from roughly 60,000 per day to about 25 per day, with the reactor 
off. 
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Figure 6: All 3 He detector signals, above the 25-count hardware threshold, occuring in the 
782 /zsec (three neutron capture times) following the event trigger, a) and b) are before any 
software cuts were applied; c) and d) after all cuts were applied. 
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Figure 7: Capture-time spectrum of neutrons detected by the 3 He counters from a simulated 
252 Cf source in the center of the target detector. The mean capture time is 265±3 /xsecs. 
Results are from 100,000 generated neutrons. 

3.2 Monte Carlo calculations 

GEANT with the GCALOR interface was used for all Monte Carlo simulations. The gcalor 
interface handles neutron transport from 20 MeV down to thermal energies. GEANT handles 
the transport of all other particles. One comparison of the data and the Monte Carlo is 
given in Figs. [?| and |[ The former shows the capture-time spectrum of neutrons detected 
by the 3 He counters from a simulated 252 Cf source at the center of the target detector. The 
mean capture time is 265 ± 3 /isecs. Fig. [8] is for the same configuration, but for real data. 
The mean capture time is 267 ± 4 /xsecs. 
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Runs 385 + 415 + 444, 2nd neutron captures 
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Figure 8: Capture-time spectrum of the second neutron detected by the 3 He counters from a 
real 252 Cf source in the center of the target detector. The mean capture time is 267±4 /zsecs. 
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3.3 Neutron Detection Efficiency 

Special neutron-calibration runs were periodically made with a 252 Cf source in the center of 
the target detector. Data from these runs were processed thru the same programs used to 
analyze the neutrino events. In particular, the same target cuts (as described above) were 
used. 

The resulting pulse-height spectra from the 3 He tubes were histogrammed for each cali- 
bration run, and the peaks fitted to Gaussians. Only those pulses within 2 standard devia- 
tions of the peak value are finally accepted as neutrons. The numbers of events with 1, 2, 3, 
4, and 5 neutrons within a given time window were tallied. The time window chosen was 3 
neutron capture times. 

Based on the known neutron multiplicity from 252 Cf fissions, one can calculate the neutron 
detection efficiency by assuming various efficiencies and comparing the observed number 
distribution with the calculated distribution. Our procedure took into account: 

• The neutron-number distribution from 252 Cf fission. 

• The neutron acceptance time window of 3 capture times. 

• The probability of an "extra" fission from the Cf source during the acceptance time 
window, which is a function of the source activity. 

This procedure yielded a mean efficiency of 0.41±0.01 for a neutron source at the center 
of the target. This value agreed well with the value derived from the Monte Carlo. As a 
result we were able to use the Monte Carlo value of 0.29±0.01 as the mean efficiency for 
single neutrons generated isotropically throughout the D 2 of the target volume. 

The efficiency for two neutrons is the square of the single-neutron efficiency (0.084±0.006). 
And the efficiency for seeing only 1 neutron, if 2 were produced is 2 x 0.29 x (1.0 — 0.29) = 
0.41 ±0.01. 

3.4 Energy calibration of inner ant is 

Since we desired to base the inner anti cut criteria on energy, both Tanks 2 and 4 must be 
energy calibrated. Periodic runs were made over the course of the experiment with a 60 Co 
source placed at various known positions in Tank 2, and beneath the center of Tank 4. (Tank 
4 was also calibrated with a 252 Cf source in that same position.) The data were compared 
with Monte Carlo simulations. Several algorithms were tested to find the best estimates of 
the energy. The best measures found were: Tank 4, sum the signals from the two PMT 
strings; Tank 2, sum the signals from the 4 vertical corner strings. Results are shown in 
Table g. 

4 Results 
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Tank 2 


Tank 4 


Uncertainty in peak position 


±15% 


±10% 


Standard deviation 


0.25E5 


0.25E3 


Max. deviation of peak over entire run 


±20% 


±10% 



Table 2: Calibration results for inner anti using a Co source. 



Reactor Up Reactor Down 



Raw 1-neutron rate (events/day) 44.62±0.59 25.28±0.68 
Raw 2-neutron rate (events/day) 2.69±0.14 1.45±0.16 

Software efficiency 0.471±0.003 0.444±0.005 

Corrected 1-neutron rate 94.66±1.24 57.00±1.53 
Corrected 2-neutron rate 5.71±0.31 3.26±0.36 



Table 3: Reactor up and down event rates. 

4.1 Event rates 

The 1- and 2-neutron event rates for both the reactor up and down data are given in Table [| 
Subtracting the reactor down rates from the up rates yields the data shown in Table [|, where 
we have also given the corresponding neutron detection efficiencies. 

The 2-neutron rate (per day) is (2.45 ± 0.48)/(0.084 ± 0.006) = 29.2 ± 6.1. To get the 
CCD rate from this value we need only correct for the effect of a nearby reactor, Reactor 
#4. It is located about 80 m from our detector. While taking data with Reactor #5 up, the 
mean power of reactor $4 was 1925 M W; while #5 was down, it was 2246 MW. This gives 
a correction factor of +0.6% to our final rates. Thus the CCD daily rate is 

R CCD = (29.2 ± 6.1) x (1.006) = 29.4 ± 6.1 

To get the NCD rate from the 1-neutron rate, two corrections must first be applied to 
the 1-neutron rate. 

• The number of CCD reactions in which only 1, instead of 2, neutrons was observed 
must be subtracted. This number is the CCD rate times the efficiency of seeing only 



1 neutron 2 neutron 
Up minus down rate 37.7±2.0 2.45±0.48 
Neutron efficiency 0.29±0.01 0.084±0.006 



Table 4: 1- and 2- neutron event rates and detection efficiencies. 
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Parameter 


Value 


% 


Detector-reactor distance 


18.5±0.1 m 


1.1 


Mass of D 2 


267.0±2.0 kg 


0.8 


No. MeV per fission 


205.0±0.7 


0.3 


Total systematic uncertainty 


1.4 



Table 5: Parameters which have significant contributions to the systematic uncertainties in 
the data rates. The last column shows the contribution of each parameter to the systematic 
uncertainty in the final event rates. 

one out of the two neutrons: 

(29.2 ± 6.1) x (0.41 ± 0.01) = 12.0 ± 2.5 

• The number of inverse-beta decays in the inner detector that leak into the target 
volume and create a single neutron must also be subtracted. From the Monte Carlo 
we estimate 22.0±0.5 inverse-beta events per day enter the target volume. Also from 
the Monte Carlo we estimate that only 5±1% of those events survive the 0.8 MeV 
inner-anti cut. Thus the number of events to be subtracted from the 1-neutron rate is: 

(22.0 ± 0.5) x (0.05 ± 0.01) x (0.29 ± 0.01) = 0.3 ± 0.1 

The corrected 1-neutron event rate is then 

(37.7 ± 2.0) - (12.0 ± 2.5) - (0.3 ± 0.1) = 25.4 ± 3.2 

Applying the single-neutron detection efficiency correction and the Reactor #4 correction 
from above, yields the daily NCD rate: 

R NCD = (25.4 ± 3.2) x (1.006)/(0.29 ± 0.01) = 88.1 ± 11.1 

4.2 Systematic uncertainties 

The significant systematic uncertainties are given in Table [|. Other possible sources of 
systematic effects which were considered, but found to be insignificant were: the calculated 
neutrino energy spectrum and the energy-calibration effects on data cuts. 

4.3 Theoretically-expected event rates 

The rates (events per day) are given by: 

N r - 

R= A^J Nv(E„)a(E v )dE v (1) 
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Isotope 


MeV/fission 


ly J1SS 


235 u 


201.7±0.6 


4.17xl0 24 


238 tj 


205.0±0.9 


5.61 xlO 23 


239p u 


210.0±0.9 


2.10xl0 24 


241 Pu 


212.4±1.0 


3.46 xlO 23 



Table 6: The energy per fission and the mean number of fissions per day for each isotope. 

where E u is the neutrino energy, N U (E U ) the daily average neutrino energy spectrum per 
MeV, N D the total number of deuterons in the target, a(E u ) the cross section for the process, 
and r is the distance from the reactor to the detector. 

The mean neutrino energy spectrum was determined from the reactor power and the 
core "burn up," i.e. the isotopic composition of the fuel, as a function of time. The reactor 
power was obtained from reactor monitoring devices several times per day. The isotopic 
composition of the fuel rods was given to us at the beginning and ending of each reactor 
cycle of about 11 months. 

The only four nuclei of importance are: 235 U||, 238 U||, 239 Pu0, and 241 Pu@. Combining 
the data in those references with the reactor power as a function of time, both the neutrino 
energy spectrum and the conversion factor from MW-hours to total number of neutrinos was 
calculated for each day. 

The energy per fission and the mean number of fissions per day are given in Table ^| for 
each isotope. 

The data- collect ion MW-hours was calculated for every day by combining the data col- 
lection times with the reactor power level at that time. The number of deuterons was 
1.605xl0 28 . 

Combining all these factors and dividing by the number of live days yields the mean 
neutrino spectrum (neutrinos/Me V/day) as shown in Table [T[ 

There has been considerable work done on the CCD and NCD cross sections in the past 
few years. Kubodera and Nozawa review the field in Ref. 0. In their Table 1, they give the 
cross sections for both the CCD and NCD reactions from threshold to 170 MeV. They state 
that the uncertainties in the values are 5%. Using the data of Ref. |§ with Eqn. [I], yields 
R NCD = 87.2 ± 4.4 and R CC d = 30.4 ± 1.5. 

4.4 Experimental cross sections 

The average cross section per neutrino is given by 

_ f N u (E u )a(E u )dE„ 
$N v {E v )dE v 

where the integrals go from the threshold for the reaction to infinity. Combining this with 
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Energy 


N y 


Energy 


N„ 


2.0 


3.56 xlO 25 


6.0 


8.98xl0 23 


2.25 


2.99xl0 25 


6.25 


6.49 xlO 23 


2.5 


2.47xl0 25 


6.5 


4.84xl0 23 


2.75 


2.09xl0 25 


6.75 


3.55xl0 23 


3.0 


1.75xl0 25 


7.0 


2.47xl0 23 


3.25 


1.45xl0 25 


7.25 


1.58xl0 23 


3.5 


1.18xl0 25 


7.5 


l.OlxlO 23 


3.75 


9.47xl0 24 


7.75 


6.13xl0 22 


4.0 


7.54xl0 24 


8.0 


3.27xl0 22 


4.25 


5.93xl0 24 


8.25 


1.35xl0 22 


4.5 


4.52xl0 24 


8.5 


8.11xl0 21 


4.75 


3.44xl0 24 


8.75 


4.89xl0 21 


5.0 


2.68xl0 24 


9.0 


2.77xl0 21 


5.25 


2.07xl0 24 


9.25 


1.65xl0 21 


5.5 


1.57xl0 24 


9.5 


1.17xl0 21 


5.75 


1.21 xlO 24 







Table 7: Time-averaged number of neutrinos per day per MeV. Energies are at lower bin 
edge. 

Eqn. p], we get 

4nr 2 R 
° ~ Nd J Nu{E u )dE p ' 
The values obtained for the NCD and CCD cross sections are given in Table || 

4.5 Improved NCD cross section 

As stated above, the CCD events create a significant background for the NCD events, and 
this background must be subtracted. The large uncertainty in our measured CCD rate 
makes a significant contribution to the uncertainty in the NCD rate. However, we note 
that our experimentally determined rates and the theoretically-expected rates agree within 
one standard deviation of each the experimental and theoretical uncertainties. Given this 
excellent agreement, we feel that an improved value for the NCD cross section may be 
calculated by using the theoretically-expected CCD daily rate (30.4 ± 1.52) rather than 
our observed rate (29.4 ± 6.1). Repeating the procedure described above, this yields an 
improved NCD rate of 86.7 ± 7.9, and a corresponding cross section of 5.98 ± 0.54 x 10~ 45 
cm 2 per neutrino. 
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NCD 


CCD 


Rate (events per day) 


88.1 ± 11.1 ± 1.2 


29.4 ± 6.1 ± 0.4 


Reaction threshold (MeV) 


2.23 


4.03 


Neutrinos (per day) 


3.86 xlO 25 


8.01 xlO 24 


Average cross section (10 -45 cm 2 per i7 e ) 


6.08 ± 0.77 


9.83 ± 2.04 



Table 8: The rates (with statistical and systematic uncertainties), reaction thresholds, total 
numbers of neutrinos above threshold, and cross sections for the NCD and CCD reactions 
as measured in this experiment 



4.6 Calculation of (3 2 

The value of (3 2 is given by the ratio of the measured neutral current rate to the theoretically 
expected rate. Thus we find 

- 88.1 ±11.1 ±1.2 

P 2 = — = i.oi ± 0.16 

H 87.2 ±4.4 

Using the improved NCD cross section determined above, we get an improved f3 2 of 0.99 
± 0.10. 

4.7 Neutrino oscillations 

Another aspect of this experiment is its ability to explore neutrino oscillation by measuring 
the ratio of the CCD to NCD rates. At reactor neutrino energies, there is insufficient energy 
to create leptons more massive than the electron. Therefore, if neutrino oscillation occurs 
at a significant level, a deficit of charged-current events compared to neutral-current events 
should be seen. This leads us to define the ratio R, where 



n ~ CCD th ' ^ ' 

NCD th 

a ratio of ratios of experimentally determined reaction rates to theoretically expected reaction 
rates. A deficit of charged current reactions could imply that some electron antineutrinos 
have oscillated to a different flavor or helicity state, either of which would imply new physics. 
We find 

' 2 ' xl±, ' A 0.334 ±0.080 



R = 88 ^nl L1 = i:z::i — = 0.96 ± 0.23 
§£§ 0.348 ± 0.004 

The error of 1% in the theoretical ratio is taken from reference ||. The neutrino-oscillation 
exclusion plot resulting from this value of R is shown in Fig. ||. 
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Figure 9: The neutrino-oscillation exclusion plot corresponding to our value of R, the ratio 
of the observed to expected ratios of the CCD to NCD rates. The solid line is the 90% 
confidence level contour; the dashed, the 95%. 
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4.8 Possible extension of this technique 

Since the theoretical error in the ratio is quite small, a high statistics, good precision mea- 
surement of R should be possible. This measurement has the potential of reaching small 
values of sin 2 29. 

In the current experiment, the CCD measurement is handicapped by the requirement 
that we observe two neutrons. The efficiency for observing this goes as the square of the 
single-neutron detection efficiency and so is necessarily small. Another method, which we 
explored but did not pursue, employs the addition of a small amount (approximately 10%) 
of light water into the heavy water target. This small addition does not effect the neutron 
detection efficiency appreciably and gives one the opportunity to observe the charged current 
reaction on the proton. Since the CCP reaction has a much larger cross-section than the 
CCD reaction, a threshold of 1.8 MeV, closer to the CCD threshold and since it can be 
detected by searching for a single neutron, one can determine the ratio of NCD to CCP with 
higher precision. 

5 Discussion 

This experiment was an improved version of our experiment done at Savannah River in the 
late 1970s. The primary improvements were in the cosmic-ray shielding, which cut that 
background by a factor of six, and an improved data-collection system. 

During the past 20 years great progress has been made in calculating the CCD and NCD 
cross sections, and they agree well with the results of this experiment. 
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